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ABSTRACT 

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way 
(MW) to be fM3i,tan < 34.4 km s _1 and the radial velocity to be in the range VM3i,rad = —109 ± 
4.4km s (Ivan der Marel et al.ll2012l ). In this study we use a large volume high resolution N-body 
cosmological simulation (Bolshoi) together with three constrained simulations to statistically study 
this kinematics in the context of the ACDM. The comparison of the ensembles of simulated pairs with 
the observed LG at the 1-er level in the uncertainties has been done with respect to the radial and 
tangential velocities, the reduced orbital energy (etot), angular momentum (fcrb) and the dimcnsionless 
spin parameter, A. Our main results are: (i) the preferred radial and tangential velocities for pairs in 
ACDM are v r = — 80± 20km s _1 , v t = 50± 10km s _1 , (ii) pairs around that region are 3 to 13 times 
more common than pairs within the observational values, (iii) 15% to 24% of LG-like pairs in ACDM 
have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the 
same sample show similar values in the inferred dimcnsionless spin parameter. It follows that within 
current observational uncertainties the quasi-conserved quantities that characterize the orbit of the 
LG, i.e. etot, 'orb and A, do not challenge the standard ACDM model, but the model is in tension 
with regard to the actual values of the radial and tangential velocities. This might hint to a problem 
of the ACDM model to reproduce the observed LG. 

Subject headings: galaxies: kinematics and dynamics, Local Group, methodsmumerical 



1. INTRODUCTION 

The Milky Way (MW) and Andromeda galaxy 
(M31) are the dominant galaxies in the Local Group 
(LG). Astronomical observations of their mass dis- 
tribution impose constraints on the standard cosmo- 
logical model . The satellite overabun dance problem 

tidal disrup- 
and the disk 



(jKlypin et all Il999ri iMoore et al.| 
tion features (McConnachie et al. 



1999) 



2009) 



dominated morphology (|Kazantzidis et al.|[2008T) are ex- 
amples of how LG studies are linked to the cos- 
mological context. Detailed studies on the Magel- 
lanic Clouds dynamics and their possible link to M31 
add to the interest of understandi ng the details of 
the LG kinematics and dynamics dBesla et ahl 120071: 
Tollerud et all 12011): iKnebe et all 120111: iFouauet et al.l 



20121 : iTevssier et al. I I2012D . However, a general concern 



in the use of the LG as a too l for n e ar-field cosmology 
dFreeman fc Bland-Hawthorn! 120021 : iPeebles fc Nusserl 
120101 ) is how typi cal is the LG reg a rding the prop- 



erties of interest (iBusha et al.l 12011 



iLiu et all l20Til : 



iForero-Romero et al.ll201lUPurcell fc Zentnerll2012D~ 

A new valuable piece of information in this issue is the 
recent observational determination of the proper-motion 
measurements of M31, which until recently had been out 
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of reach (jvan der Marel et al.ll2012h . The reported mea- 
surements set an upper bound for the tangential velocity 
of M31 with respect to the MW of tWM3i < 34.4 km 
s -1 . Together with the values of the relative radial ve- 
locity of t>rad,M3i = —109 ±4.4 km s _1 observations show 
that the relative motion of the MW and Andromeda is 
consistent with a head-on collision. With this informa- 
tion it is possible to quantify how common such a kine- 
matic configuration is in a ACDM Universe. 

This Letter presents such a study. We use a large vol- 
ume, high resolution dark matter only N-body simula- 
tion in the concordance ACDM cosmology to find a set 
of halo pairs with similar characteristics as inferred in the 
LG. We quantify these results in terms of the number of 
pairs with given radial and tangential velocities in the 
galactocentric rest frame. We also find the pairs that are 
consistent with a head on collision in terms of the ratio 
of the radial to tangential velocity ft = Vtan/^rad < 0.32 
and present these results in terms of the reduced angular 
momentum, mechanical energy and dimensionless spin 
parameter. 

In addition we make use of three constrained N-body 
simulations which are constructed to reproduce the ob- 
served large scale structure of the Local Universe on 
scales of a few tens of Mpc. The special feature of these 
simulations is that each volume features a pair of halos 
with the right characteristics to be considered LG-like 
objects. 

This Letter is structured as follows. In the next section 
we present the N-body simulations and the criteria we 
use to select LG-like halo pairs. In Section [3] we present 
the results for the dynamics in these pairs in terms of 
the tangential/radial velocities and the orbital angular 
momentum/mechanical energy. In the same section we 
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summarize these dynamical results in terms of the dimen- 
sionlcss spin parameter of the pairs. Finally, in the last 
section we comment and conclude about the implications 
of these results in the context of the ACDM model. 

2. SIMULATION AND PAIR SAMPLES 
2.1. The Bolshoi and Constrained Simulations 

The Bolshoi simulation follows the non-linear evolu- 
tion of the dark matter density field using N-body tech- 
niques. The simulation has a cubic comoving volume of 
(250/i _1 Mpc) 3 , sampled with 2048 3 particles. The cos- 
mological parameters used in the simulation are Q m = 
0.27, tt A = 0.73, a s = 0.82, h = 0.70 and n = 0.95, cor- 
responding to the matter density, vacuum energy den- 
sity, the normalization of the power spectrum, the di- 
mensionless Hubble constant and the index of the slope 
in the initial power spectrum. This set of parameters 
is compatible with the analysis of the seventh year of 
data from the Wilkinson Mic rowave Anisotropy Probe 
(WMAP) (jJarosik et al.H2011[ ). A detailed de s cripti on of 
this simulation can be found in lKlvpin etafl (|2011l ). 

With these parameters the mass per particle is m p = 
1.4 x 1O 8 /i _1 M0. In this paper we use halos obtained 
through the Bound Density Maxima (BDM) algorithm 
(jKlvpin et al.lll999al ). The halos are selected to have an 
overdensity of 200 times the critical density. Further- 
more, we only include in the analysis halos whose cen- 
ter is located outside the virial radius of any other halo. 
We have obtained the data through the public available 
Multidark database (|Riebe et al.ll2011[) . The database 
allows us to obtain the comoving positions, peculiar ve- 
locities and masses for all the halos in the simulation 
volume at z = 0. The positions and velocities of these 
haloes correspond to the average values of the 250 most 
bound particles. The Hubble flow is taken into account 
to convert the peculiar velocities into physical velocities 
and allow for a comparison with observations. We have 
verified that the main conclusions of this paper hold in 
the case of halos defined by a FOF algorithm with a link- 
ing length 0.17 times the mean interparticle distance. 

The constrained simulations we use in this Letter 
are part of the Constrained Local UniversE Simulations 
(CLUES) Q project whose main objective is to reproduce 
the large scale structure in the Local Universe as accu- 
rately as possible. The algorithm and observational con- 
str aints to construc t the i nitial conditions are described 
in iGottlober et al.l (|2010f ). We use three dark matter 
only simulations, each has a cubic volume of 64/i _1 Mpc 
on a side, with the density field sampled with 1024 3 par- 
ticles. The cosmological density parameter is J7 m = 0.28, 
the cosmological constant S1a = 0.72, the dimensionless 
Hubble parameter h = 0.73, the spectral index of the 
primordial density perturbations n — 0.96 and the power 
spectrum normalization ag = 0.817, also consistent with 
WMAP 7th year data. 

2.2. Two samples of LG-like pairs 

Based on the BDM catalogs in the Bolshoi simulation 
we construct a halo pair sample with the dynamical prop- 
erties consistent with those of the MW and M31. The 

6 http : //www . multidark . org/Mult iDark/ 

7 http://www.clues-project.org 



criteria we impose to define a LG-like halo pair are the 
following: 

1. Each halo has a mass in the range 7 x 10 11 M Q < 
M h < 7 x 1O 12 M . 

2. With respect to each halo, there cannot be any 
other halo within the mass range 7 x 10 11 M Q < 
Mh < 7 x 10 12 M Q closer than its partner. It means 
that there cannot be ambiguity on the identity of 
the pair members. 

3. The relativ e radial velocity between t he two halos 
is negative (jvan der Marel et al.ll2012T ). 

4. The distance between the center of mass of the ha- 
los must be less than l.OMpc dRibas et al] 120051 : 
van der Marel fc Guhathakurtall2008ft . 

5. There cannot be halos more massive than 7 x 
lO 12 M within a radiu s of 3Mpc with respect 
to every object centre dKarachentsev et al.l 120041 : 
iTikhonov fc Klvpinll2009D . 

6. There cannot be halos more massive than 7 x 
10 13 M Q within a ra dius of 4Mpc with respec t to 
every object centre (jKarachentsev et al~ll2004h . 

Throughout this Letter we refer to this sample as the 
full sample. This sample in the Bolshoi simulation has 
1923 pairs. Additionally, there is a sample of three (3) 
pairs constructed from the three constrained realizations. 
These pairs fulfill all the above mentioned conditions and 
additionally are located in a place with the right dis- 
tances with respect to the Virgo cluster in the simula- 
tion. 

The full observational characteristics that we take in 
this Letter for the MW-M31 pair are listed in Table [TJ 
A more reduced sample from the full sample has been 
constructed so as to obey the observational bounds on 
the masses and separation of the two main halos. These 
amount to: 

1. The separation between the center of mass of the 
halos is in the range 700 — 800kpc dRibas et all 
[20051 : Ivan der Marel fc Guhathakurtall2008f l~ 

2. The total mass of the two halos is in th e range 
1 - 4 x lO 12 M (|van der Marel et al.ll2012fl . 

Including these conditions the full sample is reduced 
from 1923 to 158 pairs. We refer to this sample as the 
reduced sample. Note that only one constrained LG-like 
object is included in this subset. 

3. RESULTS 

3.1. Radial and tangential velocities 

Figure [1] summarizes the central finding of this Letter. 
Most of the pairs in the two samples constructed from the 
Bolshoi Simulation have radial and tangential velocities 
notably different from the observational constraints at 
the 1-cr level. 

The most probable radial and tangential velocities in 
ACDM are summarized in Table [2j For the full sample 
we have «r a d,ACDM = -70 ± 10 km s _1 and n t an,ACDM = 
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Fig. 1. — Upper panels: 2D histograms of the radial and tangential velocities for LG-likc halo pairs in the Bolshoi simulation. Lower 
panels: 2D histograms of the orbital angular momentum (£ or b) and mechanical energy (etot) per unit of reduced mass calculated considering 
the halos as point masses. Left (right) panels correspond to the full (reduced) sample. The bar indicates the number of pairs in each cell. 
The loeii for the high density regions in the v T — vt plane are listed in Table [3] The half ellipse in the v T — vt plane corresponds to the 1-a 
uncertainties in observations. Equivalently, the contour line in the etot — 'orb plane encloses 68% of the Monte Carlo generated points from 
the observational values summarized in Table 1. The circles represent the positions of the pairs from the constrained simulations. 



50 ± 20 km s . For the reduced sample i> ra d,ACDM = 
— 90 ± 10 km s _1 and ftan,ACDM = 50± 20 km s _1 , where 
the uncertainties in these values reflect the minimum grid 
size needed to obtain robust statistics for the 2D his- 
togram. 

The number of pairs compatible with LG observations 
at the 1-a level are listed in Table [3] (columns 2 and 
3). In the same table (columns 4 and 5) we summarize 
the number of pairs around the ACDM values within 
the same range of absolute observational uncertainty (i.e. 
Ctan = 17km s _1 and er ra d = 4km s _1 .) 

From these results we infer that the pairs around the 
preferred phase space for ACDM are at least 13 times 
more common than pairs with the observed velocities for 
the LG. We highlight that this is a lower bound given 
that in the reduced sample none of the pairs is found in 
the interval allowed by observations. The high eccentric 
orbit of the observed LG constitutes an unlikely config- 



uration for the ACDM LG-like objects. This holds for 
the full and the reduced sample of objects. These con- 
clusions are valid at 1-cr level in the observational uncer- 
tainties. An equivalence in the abundance between the 
pairs around the observational values and those around 
the preferred ACDM velocities is reached only at the 4-cr 
level on the observational uncertainties. 

From Figure [1] it is also clear that there is a signif- 
icant number of pairs with a high tangential-to-radial 
velocity ratio. The peak in the pair number density is lo- 
cated around a region of f t = w tan ACDM/«radACDM ~ 0.7, 
while the observations suggest / t < 0.32. As summa- 
rized in Table [3] we find that only between 8% to 12% 
of the pairs are consistent with the observational con- 
straint. The three pairs from the constrained realiza- 
tions are also show higher / t ratios than observations: 
ft = 0.35,0.45,0.73. 
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TABLE 1 

Summary of kinematic observational constraints. 



^M31 ,rad 


( km s" 1 ) 


_ 1 nq 3 + 4 4 


fM31,tan 


(km s- 1 ) 


< 34.4 


^M31 


(kpc) 


770 ± 40 


r M31 


(kpc) 


(-378.9,612.7, -283.1) 


fr,M31 


(kpc) 


(-18.9,30.6,14.5) 


VM31 


( km s 1 ) 


(66.1, -76.3,45.1) 


°v,M31 


( km s -1 ) 


(26.7, 19.0,26.5) 


^200, MW 


(10 12 M Q ) 


1.6 ±0.5 


A^200,M31 


(1O 12 M ) 


1.6 ±0.5 


A^200,MW + A^200,M31 


(1O 12 M ) 


3.14 ±0.58 


log 10 A 




-1.72 ± 0.07 



Notes: 

1. The kinematic prop erties for M31 are reported in the galac- 
toccntric rcstframc ( van dor Marel et al. 2012). 

2. Values in parenthesis correspond to vector components. <T X 
represents the uncertainty on the components of vector x. 
The uncertainties correspond to l-cr values. 

3. The values for the individual halo masse s are c onsistent with 
the priors used by van dcr Marcl ct al. (2012). 

4. The observational uncertainties in the position vector cor- 
respond to a 5% in each component consistent with 
the 1-a uncertainties in the distan ce (see references in 
van der Marel fc Gu hathakurta 200|). 

5. The value for log 10 A is obtained in this Letter from a Monte 
Carlo simulation as described in Section 3. 



TABLE 2 

Summary of results from the ACDM Bolshoi simulation 



Full Sample Reduced Sample 

fM3i,rad (kms" 1 ) -70 ± 10 -90 ±10 

U M 3l!tan (kms- 1 ) 50 ± 10 50 ± 10 

log 10 A -1.47 ±0.13 -1.34 ±0.12 



Note: the velocity uncertainties correspond to the minimum bin 
size required to obtain robust statistics in Figure [T] 



3.2. Reduced Angular Momentum and Energy 

The two-body problem of point-like masses can serve 
as a proxy for the dynamics of the LG and can be used 
as a tool for studying the LG within the framework of 
the standard cosmological model. Within the model the 
dynamics of the LG is governed by the (center of mass) 
reduced angular momentum / or b = |rM3i X VM3i| and 
reduced energy e to t = 5VM31I 2 - GM/\r U 3i\, where the 
total mass is M ~ mM3i + ihmw; the reduced mass is 
/.t = mM3imMw/M and G is the gravitational constant. 

This formulation has a clear theoretical advantage if 
one considers the angular momentum and the mechani- 
cal energy as quasi-conserved dynamical quantities. This 
means that, after some formation time, these quantities 
do not significantly evolve. 

However, this formulation has an observational disad- 
vantage. The reduced energy and angular momentum 
are derived from very different kinds of observations, in- 
creasing the uncertainties in their final determination. 

We use Monte Carlo sampling to estimate the reduced 
energy and angular momentum from the observed prop- 
erties in the LG listed in TablcJTJ Lower panels in Figure 



Q] also presents the results in the plane e to t — ^ortfl As 
expected, the e to t and l OI ^ constraints are less restrictive 
compared with the radial and tangential velocity con- 
straints; 15% (24%) of pairs in the full (reduced) sample 
obey the 1-a observational constraints. We note that for 
the reduced sample the preferred region is outside the 
1-a observational contours. 

All the constrained LGs of the full and reduced sam- 
ple are consistent with the 1-a observational uncertainty. 
However a an hypothetical increase by a factor of 2 in ac- 
curacy of the tangential velocity, showing that it is below 
17 km s" 1 , would bring these pairs and the ACDM ex- 
pectation outside this uncertainty region. 

3.3. Dimensionless Spin Parameter 

The dimensionless spin parameter A (|Peebleslll971f ) is 
used here to characterize the dynamical state of the ob- 
served and simulated LGs. The parameter measures the 
dynamical role of the angular momentum in terms of the 
gravitational attraction and is defined by 

_ /J 3/2 ^orb\/etoT , , 

We compare the distribution for A obtained from the 
pair population in the Bolshoi simulation and the dis- 
tribution from the Monte Carlo simulation used to esti- 
mate the uncertainties on etot and Z or h. Figure [5] shows 
the likelihood distribution of the observational errors and 
the spin parameter of the LG-likc pairs of the full (left 
panel) and the reduced (right panel) samples. The value 
of A of the constrained LGs is shown as well. 

The expected observational value for A is slightly in- 
consistent with the statistics derived from halo pairs in 
ACDM. The spin parameter of the simulated LG-likc 
objects is skewed towards higher values compared with 
the observational estimation. Inspection of Figure [2] also 
shows that the A distribution is close to a log-normal one 
and thus resembles the distr ibution of the sp in parameter 
of DM halos in simulation dBett et al.ll2007D . 

The estimated value of log 10 A of the observed LG is 
— 1.72 ± 0.07, lower than the mean value of log 10 A of 
-1.47 ±0.13 (-1.34 ±0.12) for the full (reduced) sample. 
The values of the spin parameter of the 3 constrained 
LGs are log 10 (AcLUEs) = —2.21, —1.72 and —1.65, with 
the latter value belonging to the sole constrained LG of 
the reduced sample. 

4. CONCLUSIONS 

We have presented a comparison between the observed 
kinematics for the M31 in the galactocentric rest-frame 
and the expectations for a large N-body cosmological 
simulation in the ACDM cosmology. In the simulation 
we select a large sample of pairs that fulfill isolation crite- 
ria for the LG. We select a sub-sample which obeys more 
stringent observational constraints on the mass and sep- 
aration of the MW and M31. 

While the observations show that M31 moves towards 
us on a highly eccentric orbit, the simulation shows that 
the most common configuration at z = has values 

8 The l-cr contour in the lower panels of Figure [l] docs not cross 
the zero level, as can be naively expected from Figure [T] This 
happens for the si 80% contour level. 
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TABLE 3 

Summary of the comparison of the observational results against ACDM. 



Physical 


(%) Pairs consistent 


(%) Pairs consistent 


(%) Pairs with highest 


(%) with highest 


property 


with observations (1-cr) 


with observations (1-ct) 


likelihood in ACDM 


likelihood in ACDM 




(full sample) 


(reduced sample) 


(full sample) 


(reduced sample) 


Vi-Vt 


(0.4%) 8/1923 


(<0.6%)0/158 


(1%)23/1923 


(8%)13/158 


Ctot-^orb 


(15%)298/1923 


(24%)38/158 






logi A 


(13%)257/1923 


(9%)15/158 






rt = vt/vt 


(12%)242/1923 


(8%) 13/158 







Notes: 

1. In columns 2 and 3 consistency is defined from 1-cr uncertainties for v T -Vt, et t-£ rb and log 1( ) A. In the case of rt, is defined from 
the constrain rt < 0.32 derived from the I— a observational uncertainties in the radial and tangential velocities. 

2. In columns 4 and 5 the number of pairs around the preferred region in v T — vt in ACDM (as shown in Fig. [T] and summarized in 
Table [2]l are calculated using the same observational 1-cr absolute uncertainty on v r and vt- 




Fig. 2. — Normalized histograms of Peebles' spin parameter A for the pairs in the Bolshoi simulation and its inferred values for the LG 
from the observational constraints from a Monte Carlo simulation. The left (right) panel corresponds to the full (reduced) sample. The 
vertical lines with the white dots represent the values inferred for pairs in the constrained simulations. 



Wrad,ACDM = -80 ± 20 km S 1 and W t an,ACDM = 50 ± 10 

km s , where the error bars are estimated to include 
the results from the two pair samples. 

Using the same absolute values for the uncertainty in 
the observed velocity components, wc find that none of 
the halos in the reduced sample are compatible with ob- 
servations. This makes pairs with the preferred ACDM 
values at least 13 times more common than pairs com- 
patible with the observational constraint. Additionally, 
pairs with a fraction of tangential to radial velocity 
ft < 0.32 (similar to observations) represent 8% to 12% 
of the pairs. 

Approximating the LG as two point masses we express 
the above mentioned results in terms of the orbital an- 
gular momentum / or h and the mechanical energy etot per 
unit of reduced mass. We find that the uncertainties in 
the tangential velocity, the square of the norm of the ve- 
locity and the total mass in the LG are less constraining 
on the number of simulated pairs that arc consistent with 
the observations. Nevertheless, in the case of the reduced 
sample there is a slight tension between simulation and 
observation. A reduction by a factor of 2 in the obser- 
vational uncertainty on the radial velocity would clarify 
this issue. 

We also use the A spin parameter to gauge the dynam- 
ical state of pairs. The values for the orbital angular 



momentum and energy, merged into the A spin parame- 
ter, are in mild disagreement with the observational con- 
straints. 

In the three pairs from constrained simulations we find 
kinematics dominated by radial velocities. However their 
velocity components differ from the observational con- 
straints and their mechanical energy and orbital angular 
momentum are in broad concordance with observations. 
There is only one pair that fulfills all the separation, to- 
tal mass constraints and matches the most probable value 
for the dimensionless spin parameter A inferred from ob- 
servations. 

Summarizing, we see a broad agreement in the total 
angular momentum and energy and a marked difference 
with the precise balance between today's radial and tan- 
gential velocities, i.e. head-on collisions of MW-like ha- 
los, as described by 1-cr uncertainties in observations, are 
not common in ACDM. 

Under the approximation of conservation of the orbital 
angular momentum and mechanical energy, this could 
only be explained if the initial conditions for the forma- 
tion of the Local Group are special in comparison to the 
initial conditions of any other pair of dark matter halos 
in the ACDM cosmology. Investigating this perspective, 
by a thorough characterization of the small pair sample 
consistent with observations in terms of halo properties 
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such as concentration , spin and place in the cosm ic web 
(jForero-Romero et al.|[2009t iHoffman et al.ll2012ft is un- 
derway (Bustamante et al. in prep.). 

The results presented in this Letter open a new window 
into the question of how unique, if at all, is the LG in 
a cosmological context. This will continue to be studied 
within the framework of the CLUES project. 
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